Coherent Umklapp Scattering of Light from Disordered Photonic Crystals 
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A theoretical study of the coherent light scattering from disordered photonic crystal is presented. 
In addition to the conventional enhancement of the reflected light intensity into the backscattering 
direction, the so called coherent backscattering (CBS), the periodic modulation of the dielectric 
function in photonic crystals gives rise to a qualitatively new effect: enhancement of the reflected 
light intensity in directions different from the backscattering direction. These additional coherent 
scattering processes, dubbed here umklapp scattering (CUS), result in peaks, which are most pro- 
nounced when the incident light beam enters the sample at an angle close to the the Bragg angle. 
Assuming that the dielectric function modulation is weak, we study the shape of the CUS peaks for 
different relative lengths of the modulation-induced Bragg attenuation compared to disorder-induced 
mean free path. We show that when the Bragg length increases, then the CBS peak assumes its 
conventional shape, whereas the CUS peak rapidly diminishes in amplitude. We also study the sup- 
pression of the CUS peak upon the departure of the incident beam from Bragg resonance: we found 
that the diminishing of the CUS intensity is accompanied by substantial broadening. In addition, 
the peak becomes asymmetric. 
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I. INTRODUCTION 

Since the first experimental observationsETH, the phe- 
nomenon of coherent backscattering (CBS) of light from 
disordered media has been the subject of intense theo- 
retical and experimental studiesEll (see also Ref. for 
the most recent review). The underlying mechanism for 
the CBS was identified as interference of clockwise and 
counterclockwise scattering paths. This was understood 
already in the early works by analogy to weak localiza- 
tion of electrons. It has been also pointed outEH that this 
picture fully captures the physics of the coherent scat- 
tering only if there are no forbidden directions for the 
propagation of light in the absence of disorder. These 
forbidden directions emerge in systems with periodic 
spatial modulation of the dielectric function, or in-other 
words, in photonic crystals with incomplete gapstL3. In 
the presence of periodicity, the enhanced scattering of 
light may occur not only in the backscattering direc- 
tion, but in other directions as well. Roughly speaking, 
the additional peaks in the scattering intensity can be 
regarded as periodicity-induced diffraction satellites of 
the CBS peak. Their origin is illustrated in Fig. 1. In 
the presence of the periodic modulation of the dielectric 
function, the light wave vector is determined only up 
to the vector cr of the reciprocal lattice. Hence, upon 
entering the medium, light with wave vector k acquires 
a satellite component with y projection of the wave vec- 
tor ^ — ' ^ 

sin 9 — a. Importantly, the same argument 
also applies for the coherently backscattered light with 
wave vector — k, when it propagates inside the medium 
on the way out, namely it also acquires a component 
with y projection k y = — fcsin# + a (Fig. 1). This com- 
ponent gives rise to a satellite of the CBS peak in the 
direction 9' where sin 9' = (a — ksin9)/k. We dub this 




Fig. 1. Schematic plot of light diffraction in a photonic 
crystal. The diffraction satellites of the incident (at angle 9) 
and outgoing plane waves are shown with dashed lines. The 
diffraction satellite of the backscattered wave (-k) results in 
outgoing wave emerging at an angle 8' . 

peak in the reflected light intensity as coherent umklapp 
scattering (CUS). We note that the above picture is only 
illustrative. In reality, in addition to the process shown 
in Fig. 1, a variety of diffraction processes contribute to 
the formation of the CUS peaks. In general, the total 
number of the CUS peaks is determined by the number 
of reciprocal lattice vectors for which — k + cr\ < k. 

It follows from the above qualitative picture that the 
magnitude of the CUS peak is governed by the ratio of 
the disorder-induced mean free path, I (the elementary 
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step of light diffusion) and the characteristic length, £, of 
formation of the diffraction component. This formation 
occurs most efficiently for the incidence angles = 8 B , 
corresponding to the Bragg condition k y = cr/2, i.e. 
s'm9 B — a /(2k). In this case £ coincides with the Bragg 
decay length, L B , which is the decay length of evanescent 
wave with frequency in the middle of the photonic stop 
band. Away from Bragg resonance, the length £ increases 
resulting in suppression of the CUS peak. We note that 
the limiting case L B — > and = B was considered in 
Rcf. [l3]. In this case the CBS and CUS peaks are simply 
the mirror images of each other. In the case of photonic 
crystal with incomplete gap, L B is large (kL B 3> I). The 
question of interest is then: how do the magnitude and 
shape of the CUS peak depend on L B /l and on the de- 
viation of the incident beam from the Bragg angle. This 
question is studied in the present paper. We generalize 
the approach of Refs. || and [li] to the case of periodi- 
cally modulated dielectric medium and derive analytical 
expressions for the CBS and CUS peaks. We trace the 
evolution of the amplitude and width of the CUS peak 
with changing the parameter L B /l and the incidence an- 
gle, 9. We found that the most pronounced variation 
occurs within the domain L B < I and (9 — 9 B ) < (kL B )~ l . 

The paper is organized as follows. In Section II we 
review the conventional derivation of the CBS and in- 
troduce the modifications necessary to take into account 
the Bragg reflection in periodic structure. In Sec. Ill 
we derive analytical expressions for the CUS and CBS 
albedo from disordered photonic crystal and analyze lim- 
iting cases. In concluding remarks, Sec. IV, we outline 
generalization of our theory to the case of arbitrary di- 
rection of the modulation wave vector. 



malized mean incident field at point r inside the medium, 
and U (r 1; r[, r 2 , r'p^is the sum of all scattering diagrams 
with ends strippedtlJ. 

There are two leading contributions to 
U(r 1 ,r' 1 ,r 2 ,r' 2 ). The first one comes from the ladder 
diagrams and describes the background incoherent scat- 
tering. The second contribution, which is responsible for 
the coherent enhancement of scattered intensity, repre- 
sents the sum of maximally crossed diagrams 

U * = ^ p (TxSx)&{ri f a )*(rj - r 2 ), (2) 

where P(r l ,r' 1 ) is the stationary probability distribution 
to travel diffusively from r to r[ inside the scattering 
medium. For the disordered medium with a boundary, 
this probability distribution is given by 



P(r v 



1 



1 



r = 



AnDlr, 



(3) 



where D — lc/3 is the light diffusion constant and c is the 
speed of light. The first term in Eq. (|^) is a conventional 
propagator in the bulk medium. The second term en- 
sures the boundary condition P(z ,r' 1 ) = P(r l7 z ) = 
that is imposed to describe the diffusion inside a semi- 
infinite medium. The point r* is the mirror image of the 
point r\ with respect to the trapping plane located at 
zq w — 0.7Z (see Ref. and references therein). 

The mean incident field in the case of translationally 
invariant (on average) system is given by 



ty iac (ri) = exp ^~ 



z\ 



I cos 9 



(4) 



II. COHERENT ALBEDO 
A. Coherent backscattering from disordered medium 

In order to introduce the notations used throughout 
this paper, we first review the conventional derivation 
of the CBS Mlowing Ref. [l4|. Neglecting the polar- 
ization effectsEj, the intensity I(R) reflected from the 
medium when illuminated with incident flux Fo and 
which is observed at point R, is given by the well known 
expression^, 



I(R) = F / G(R, r' 1 )G*(R, r' 2 )U(r 1 ,r' 1 ,r 2 , r' 2 ) 



x *mo(ri)*L( r 2) d 3 r 1 d 3 r 2 d 3 r[d 



V 3l 



(1) 



where G(R, r) is the mean propagator from r (inside the 
medium) to the observation point _R, ^ , inc (r) is the nor- 



where k is the wave vector of the incident light ( | k \ = k = 
lo/c) and 9 is the angle of incidence (Fig. 1). The first 
term in the exponent describes the decay of the incident 
mean field amplitude due to scattering. 

For the observation point R in the far field region, the 
asymptotic expansion of Green's function G(R, r±) is 



G(R, ri 



4vri? 



exp 



I cost 



ik r% 



(5) 



where fe' is the wave vector in the direction of observa- 
tion, \k I = uj/c, and 8' is the angle between fc and z 
axis. 

Substitution of Eqs. (§|)-(||) into Eq. (|f|) results in the 
following expression for the CBS albedo, a, defined as 
the scattered intensity divided by the incident flux and 
the sample area S, 



&{k, k') = ^2 S i3 J d 2 pdzidz[ exp 



iqp + in(z\ 
1 



-b 



(zi + z[) 



I 



Vp 2 + ( Zl - z[y Vp 2 + (*i + A + 2z o) 2 



(6) 
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where p is the component of (r\ — r\) paral- 
lel to the medium boundary, k = (k + k) z , 
q = {(k + k') x ,{k + k') y }, and b = (l/cos0 + l/cos0'). 
Due to the presence of the fast oscillating exponent, 
exp(iqp), the integral in Eq. is non-zero only within 
a narrow interval \9' — 0| ~ (&0 1 around the backscat- 
tering direction, i.e. q = 0. In Ref. 14 this integral was 
evaluated analytically for small angles 0, 0'. In fact, a 
general expression that is valid for arbitrary 0, 0' can be 
obtained, 



a(fc,fc') = f(n,q,b) 



1 



8ttS (k/) 2 
1 - exp(-2qz ) 



(6 + g/) s 



(7) 



In the next subsection we trace how the albedo ([?]) is 
modified in the presence of a weak periodic modulation 
of the dielectric function inside the scattering medium. 



B. Coherent scattering in the presence of a photonic 
crystal 

1. Modification of the wave amplitudes 

Our main observation is that incomplete photonic gap 
affects light diffusion only weaklyll3. The probability that 
a certain step of the diffusive motion occurs in the for- 
bidden direction can be estimated as (fcL B ) _1 <C 1. On 
the other hand, before the first and after the last scat- 
tering events the light propagates along fixed directions 
and 0'. If either or 0' is close to the forbidden di- 
rection, then light propagation will be strongly affected 
resulting in a significant change in the albedo. This ob- 
servation suggests that in order to calculate the coherent 
albedo from a photonic crystal with incomplete gap, it is 
sufficient to modify only the incident field amplitude in 
Eq. (4) and the Green function of emerging light in the 
Eq. (5), without changing the propagator in Eq. (||). 

To model the incomplete band gap, which is narrow 
compared to the Bragg frequency one can keep only a 
single harmonics in the spatial modulation of the dielec- 
tric function, 



Se(y) = 2Se cos(ay). 



(8) 



Here we assumed for simplicity that the direction of mod- 
ulation is parallel to the boundary as shown in Fig. 1. 
Generalization of the results to an arbitrary angle be- 
tween the boundary and the modulation wave vector is 
outlined in the concluding remarks. 

As was pointed out above, the Bragg resonance con- 
dition has the form k y — a/2. For a wave propagating 
in a boundless medium this condition would lead to an 
amplitude decay oc exp[— Im(k y )y] in the y direction, for 
light within the frequency range Au acde (the pho- 
tonic stop band). However, the boundary conditions en- 
force a real k y value (see Fig. 1). Consequently, instead 



of causing a finite Im(fcj,), the Bragg condition manifests 
itself in splitting of the z projection of the wave vector. 
This may be seen from the following relation between the 
components of the wave vector k inside the medium in 
the vicinity of the Bragg resonance: 



cos" 



(9) 



where Sk z = k z — k cos B and Sk y = k y — k sin B . The 
derivation of Eq. (^|) is sketched in the Appendix. With 
k y = k y = fcsin0 fixed by the boundary conditions, 
Eq. (0) yields two values of k z , namely k z — k cos B ± 57 
with 



1 



2k cos 
tan B 
2L B 



— V(fffccos0 B /3) 2 + (fc 2 fe) 5 
V(2/cL B cos0 B (3f + 1 , 



(10) 



where (3 — — B is the deviation from the Bragg an- 
gle and L B = 2 sin 2 9 B /(aSe) is the Bragg length. In the 
region z > 0, the field components with z projections 
fccos0 B + f2 and fccos0 B — f2, which comprise the waves 
k and k — cr shown in Fig. 1, are coupled to each other. 
This leads to the following modification of Eq. ([|) for 



^^(r) w exp 



I cos ( 



x {c n , 4> (zy kBr + iSnA*)e i(kB ~' T)r } - ( n ) 
where the functions Cq^(z) and Sq^(z) are defined as 



Cn,4>(z) = cos(57z) — icoscj)sm(£lz), 
Sn,(t>(z) = sin0 sin(f2,z), 

and <p is determined by the relation 

1 



sin = 



(12) 
(13) 



(14) 



v/(2fcL B cos0 B /3) 2 + 1 

The first term of the field amplitude in Eq. ([ll]) is 
the wave transmitted through the interface and traveling 
along direction close to the direction k of the incident 
wave, (k B ) z = A:cos0 B , and (fc B )a;iy = k xy = fcsin0. The 
second term is the diffracted satellite wave. As is seen in 
Eq. ([l3"|), the characteristic length scale at which the lat- 
ter wave is formed is £ = 1/S7. In the limit Lb — > oo the 
incident wave does not change upon crossing the bound- 
ary. Indeed, the satellite wave in Eq. ([ll]) vanishes due 
to the sm.(j> prefactor in Eq. (|l3|), whereas Cn^(z) turns 
into exp(— ifcsin0 B [3 z) and makes up for the difference 
between (k B ) z and k z — fccos0. 

The modification due to the Bragg scattering should 
be incorporated into the emerging wave propagator G in 
a similar way. Taking into account that after the last 
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scattering event a wave exp(ifc'r) that emerges at angle 
9' = 9 B + 59 close to the Bragg resonance also has a 
diffracted component inside the medium, Eq. (||) trans- 
forms into 



G(R,r) 



A-kR 



exp 



I cos 9 




Fig. 2. Schematic illustration of different contributions to 
the coherent albedo. CUS peak originates from (a) diffrac- 
tion of the backscattered light and (b) backscattering of the 
diffracted light. 



J 



x {Cn>^(z)e- l( - k 'B-^ r + iSn>^<(z)e- lk >} , (15) 

where SI' and <p' are given by the same expressions (|ic|), 
(|l4|), where f3 is replaced by 59. The components of 
the wave vector k' B inside the medium are given by 
(k B ) z — —kcos9 B and {k B ) X y — k' xy = ks'm9'. Expres- 
sion (|l^) is written specifically for the observation point 
in the CUS direction, i.e. (k B ) y > (Fig. 1). Similarly 
to Eq. ([Til), G(R,r) represents the sum of two terms, 
which wc dub here as C-term and S'-term, respectively. 



2. CUS and CBS albedo 



Substitution of Eqs. (O) and fll5| ) together with 
Eqs. (Q) and (^) into Eq. ([]]) produces a sum of various 
terms with fast oscillating exponential factors. The terms 
with factors e ±4<Tri or e ±lCTr i average out upon integra- 
tion. Since Eq. ( |l5| ) is valid only in the vicinity of the 
Bragg resonance k w — k + er, then nonrcsonant terms 
should be discarded. Note however, that, in addition to 
the oscillating terms, the product of field amplitudes (|ll| ) 
and propagators (|l^) in Eq. ([!]) contains now two terms 
proportional to exp[i(fe' B + k B — er)p], which do not van- 
ish if the outgoing light is parallel to the CUS direction. 
The total CUS albedo is determined by these terms and 
is thus given by 



j (fc,fc / ) = j d 2 pdz 1 dz[ exp 



iqp 



Cn,4>( z i) C Q\^i z i) s Q',<t>'( z 'i)Su ,<f>(z[)] 



(16) 



where q = (fc B + k' B - <r) xy . 

The origin of the two terms in the coherent scattering albedo is schematically illustrated in Fig. 2. In Fig. 2a the 
C-component of the incident light (solid line) first experiences coherent backscattering and then is diffracted into the 
S component (dashed line). In Fig. 2b the C-component of the incident light is first diffracted into the ^-component, 
the backscattering of which provides the second contribution to the CUS. 

The origins of the two contributions to the CBS can be traced from Eq. ([!]) in a similar way (see also Fig. 2). The 
only technical difference between the derivations of the CBS and CUS is that for the observation point R in the CBS 
direction^, e. k' m —k), cr should be added to the wave vectors in both oscillating exponents of the Green function 
in Eq. 0. 

Substitution of the modified propagator into Eq. (|l|) and selection of non-vanishing resonant terms produce the 
following expression for the CBS albedo 



i(k,k') = J d 2 pdz 1 dz[ exp 



iqp 



I 



P(ri,ri) 



(17) 
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with q = (fe B + k' B ) xy . In the limit Lb — > oo, the 
diffracted waves vanish, Sn,<j> — ► 0, so that the only con- 
tribution to the albedo that survives in this limit comes 
from the first term of Eq. (fij 



III. RESULTS AND DISCUSSION 
A. Analytical results 

1. Expressions for CBS and CUS albedo 

The additional oscillating factors C(z) and S(z) in the 
integrands (|l6|), ( p7| ) compared to Eq. (g) can be formally 
absorbed into the decrement b by adding to it imaginary 
parts of the type i(ft±ft')Z. Then all the contributions to 
the coherent albedo can be conveniently expressed with 
the help of an auxiliary function, f(n,p), defined as 



= f( K , Q, b + iP l ) + /(«, q, b ~ ipl) 



1 



AnSD l + X 2 



b-plX 1 - e- 2 ? 20 



b 2 + p 2 P 



ql 



, (18) 



where /(k, &) is the shape of the CBS cone given by 
Eq. (0); the parameters D and X are expressed through 
the arguments of the function / as follows 



D = ( K 2 -p 2 )l 2 + (b + ql) 2 , 
2pl(b + ql) 



X = 



D 



(19) 
(20) 



In Eq. ( |18| ) the wave vector q is equal to q = k cos 6 B \/3 — 
80\ for the CBS, and q = kcos0 B \3 + 60\ for the 
CUS; in the vicinity of the Bragg resonance we can set 
b = (coaO)- 1 + (cosO')- 1 « 2(cos(9 B )- 1 . Finally, the the 
two contributions to the CBS albedo in Eq. (17) take the 
form 



«CBS 



16 



/(o, n + n') + /(o, n - ft') 



(1 + COS 2 0)(1 +cos 2 ^') 



16 



(1 



cos 



-/(ft, ft' 



/(ft + ft',0) + /(ft-ft',0) 
(1 +cos 2 0')sin 2 - 



COS0 COS( 



(2) 
«CBS 



sm 



sin 



16 

/(o,ft-ft') 



/(ft',ft)- 

/(ft + ft',0)-/(ft-ft',0)] , 
/(ft + ft', 0) + /(ft - ft', 0) + /(0, ft + ft') + 

/(ft', ft) -/(ft, ft') . 



(21) 



(22) 



Analogously, the CUS albedo Eq. (|l|) can be expressed through the function / in the following way 



sm sm < 



{< 



1 + cos 6 cos 6') 



/(0, ft - ft') + /(ft - ft', 0) - /(ft, ft') - /(ft', ft) 



(1 - cos 4> cos 4>') /(0, ft + ft') + /(ft + ft', 0) - /(ft, ft') - /(ft', ft) | , 



(23) 



Expressions (21 )— (|23|) are our main results. Below we analyze two limiting cases of small and large L E 



r 



2. Limiting cases 

Let us trace how the conventional CBS cone is re- 
covered in the limit Lb — > oo. In this limit we have 

(2) 

4>, 4>' — > 0, so that a cus and a C Bs containing sin0 and/or 
sin 4>' as prefactors, vanish. Substituting ip — <f>' = 

into Eq. @ we obtain aQL = /(ft - 0', «?,&)• Tak- 
ing the limit L B — > oo in Eq. (^0|), we get (ft — ft') — » 
/c sm C7 B (3 — 60). Correspondingly, /(ft — ft', g, 6) reduces 
to Eq. |). 

Consider now the opposite limitliij, L B /^ <1, In this 
limit </>,<//—► 7r/2. It follows from Eq. (M), that with de- 
creasing L B both ft and ft' diverge, while (ft — ft') — > 0. 
As a result, as can be also seen from Eq. (Eq), all of the 
terms in Eqs. (pl|)-(p3|) that contain ft, ft^or ft + ft' as 



at least one of the arguments, vanish. Then it is straight- 
forward to check that both a CBS and a cus take the form 
/(0,0]/4 = /(0, q, b)/2. We thus recover the result of 
Ref. [l3| that in the limit of strong modulation, CBS and 
CUS cones are the mirror images of each other. Concern- 
ing the shape of the cones, it is given by the conventional 
expression (^) for the coherent albedo. Concerning the 
peak heights, they are two times less than the height of 
the CBS peak in the absence of modulation. 



B. Shapes of the CBS and CUS peaks 



Expressions (glJ)-(gSj) for the CBS and CUS albedo 
are valid only in the vicinity of the Bragg resonance 
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£/cos0 u 89 




£/cos9 R 8e 




k I cos9 B 86 

Fig. 3. CBS (thin line) and CUS (bold line) peaks at 
L B = 0.3/ normalized to the CBS peak height at L B /l — > oo 
are plotted for detuning (a) (3 = 0; (b) \/3\ = (kl cos 6 B ) _1 ; (c) 
|/? | = 2(kl cos9 B )~ 1 . The peak maxima occuring at 5$ = (3 
(89 = -[3) for the CBS (CUS) are shifted to 89 = for con- 
venience. 

6,6' « #b. Note, however, that only in this region CUS 
has an appreciable amplitude. To illustrate this we plot 
in Fig. 3 the CBS and CUS cones for different detunings, 
[3 = 6 — 6 B , of the incident beam from the Bragg angle for 
L B /l = 0.3. It is seen that CUS practically dies out at 
[3 « (2kL B cos0 B ) _1 <C 1. It is also seen that as the am- 
plitude of the CUS peak falls off, the peak also becomes 
asymmetric. The behavior of the CUS and CBS peak 
heights with the detuning (3 is summarized in Fig. 4. 

Consider now the case of exact resonance, (3 = 0. As 
was discussed above, the relation between the CBS and 
CUS peaks is governed by the dimensionless parameter 
L B /l. In Fig. 5 we show the heights of both peaks as func- 
tions of L B /l. It is seen that the amplitude of the CBS 
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Fig. 4. CBS and CUS peak heights normalized to the 
CBS peak height at L B /l — > oo plotted vs. the detuning 
(klcos9 B )(3 at L B = 0.3/. 

peak saturates already at L B /l > 2; at the same time, 
the CUS peak diminishes by an order of magnitude. 

In Fig. 6 we illustrate the evolution of the CUS cone 
at f3 = with increasing L B /l. For reference the CBS 
cone for L B /l = 0.3 is also plotted in Fig. 6 (solid line). 
We note that, while the CUS peak keeps narrowing for 
L B < I, the shape of the CBS cone remains practically 
unchanged with L B /l. 



IV. CONCLUSIONS 

In this paper we have studied coherent light scatter- 
ing from a disordered photonic crystal with incomplete 
band gaps. We have demonstrated that the crystal di- 
electric function periodicity gives rise to additional coher- 
ent albedo peaks in non-backscattering directions. These 
peaks emerge as the angle of incidence, 6, approaches the 
Bragg resonance, 6 ~ & B . For simplicity, the considera- 
tion was restricted to the case where the modulation wave 
vector, a, is parallel to the crystal boundary. In this case, 
and under the Bragg resonance condition 2fcsin(9 B = a, 
the direction of the additional (CUS) peak coincides with 
the reflection direction, since the CUS wave vector is 
equal to k' = —k + cr. The component of k! parallel 
to the boundary is equal to k' y = —ksind + cr. Then at 

6 = 6 B we have k' y = k y = a/2. It is important to note 
that this relation is specific for cr parallel to the bound- 
ary (and for the wave vector of incident light lying in yz 
plane) . 

Consider now the case where cr and the incident light 
wave vector still lie within yz plane, but cr is no longer 
parallel to the boundary, but rather makes an angle 

7 with the y axis. In this case the Bragg condition 
takes the form 2fcsin6> B = a, where 6 B = sin -1 (^) — 7. 



G 











L B /l 



L. 



Fig. 5. CBS and CUS peak heights normalized to the CBS 
peak height at Lb/I — > oo plotted vs. L B /l. 

It is easy to see that together with the CUS condi- 
tion k' = —k + cr this determines the CUS direction 
Ocvs — Ob + 27, whereas the reflection angle IS C7 r — Pb . 
This latter case can be exploited in experiments. 

Generalization of our theory to the case of finite 7 
amounts to the replacement of the functions C and S 
in Eqs. ©, @ for the CUS and CBS by 



Cq 7 = cos flz — i sin cj> sin Clz, 



cos 9 E 



cos(6» B + 2 7 ) 



1/2 



sin^sinfiz, 



(24) 
(25) 



where 



20 = 



sin # F 



i B cos(6> B + 27) 



(2£lL b cos(6> b +7)/3) 2 + 1, (26) 



is the modified splitting between the two solutions of 
Eq. (0) for k z . The modified Bragg length in Eq. (|2G|) is 



defined as L B = L B y cos# B / cos(# B + 27). The parame- 
ter tfi in Eq. (^) is still given by Eq. (|T|) with L B — > _L B . 

To derive Eqs. (J2^)— (|26|) it is convenient to perform a 
rotation of the coordinate system by the angle 7. Then 
the additional factors in (|25|), ( pf ) as compared to 7 = 
emerge due to the modification of the boundary condi- 
tions. Other steps of the derivation remain unchanged. 

It is seen from Eqs. (p5f), (26) that the prefactors 
diverge at 27 = ^ — 9 B . This divergence corresponds 
to the physical situation where the diffracted compo- 
nent of the incident light is aligned with the boundary. 
With regard to the coherent scattering, the condition 
27 = ? — 9 B manifests a crossover to a new regime. 



k I cos0 B 80 



Fig. 6. Normalized CUS albedo at exact resonance, (3 — 0, 
shown for L B = 0.31 (dashed line) and L B — I (dot-dashed 
line). Solid line: normalized CBS albedo at ft = and 
Lb = 0.3Z. 

In this regime the diffracted wave does not "fit" into the 
medium, so that the CUS peak is absent. Formally, 
for 27 > -| — 6 B Eq. (|9|) does not have real solutions 
for k z for small /3, which corresponds to opening of the 
photonic band gap in z direction. As a result, the CBS 
peak exhibits an anomaly in the vicinity of the Bragg 
conditior£3. 

We did not address in the present work the modifi- 
cations of the theory caused by a difference of refrac- 
tion indices on two sides of the interface. The impact of 
the light refraction on the CBS was intensively studied 
(see Ref. ^ and references therein). A non- trivial conse- 
quence of the refraction index mismatch is that in course 
of diffusion the light wave can "strike" the boundary at 
an angle exceeding the angle of total internal reflection. 
Basically, the CUS can be viewed as a mirror image of 
the CBS, hence the effect of the index mismatch on CUS 
is similar to that on CBS. 
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APPENDIX: 

In the vicinity of the Bragg resonance, k y ~ er/2, we 
employ the coupled wave approach, i.e. we search for the 
field inside the medium with dielectric function given by 
Eq. (||) in the form of a sum of two waves with wave vec- 
tors k and k — a and amplitudes and A^ ^., respec- 
tively. Substitution of this form into the wave equation 
results in a following system of coupled equations for the 
wave amplitudes 
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[(fccos(9 B ) 2 -k 2 z 
[(* 



k 2 Se A 



k 



cos 



2fcsin6» B (fc y - fcsin0 B )] A k + k 2 5e A k a = 
? B ) 2 - k 2 z + 2ksm6 B (k y - fcsin<9 B )] A fc cr = 0, 



(1) 



where we used the definition of the Bragg angle 
sin6> B = a I '(2k). The system ([!]) together with the as- 
sumption cos 9 « cos # B lead to Eq. (||) . 

It is convenient to express the modulation strength Se 
in terms of the Bragg decay length, L B , given by 



Lb = 



2 sin 2 B B 
aSe 



(2) 



It is easy to see from Eq. (||) that the meaning of Lb 
is the decay length for light of frequency lo — kc in the 
middle of the Bragg gap, when the medium boundary is 
perpendicular to the modulation direction y. 

If the condition given in Eq. (||) is satisfied, then there 
exists a nontrivial solution of the system (Q), 



Ai. 



.4 



k-cr cos^il 



.4 



k 



(3) 



with <j> defined by Eq. The two signs in the numer- 

ator correspond to the two solutions of Eq. (|9|) for k z . 
Upon matching the solution for the field amplitude with 
the incident wave, the resulting field can be cast in the 
form of Eq. (0). 
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